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AN  INVARIANCE  PRINCIPLE  FOR  DYNAMICAL  SYSTEMS  ON  BANACH  SPACE: 

APPLICATION  TO  THE  GENERAL  PROBLEM  OF  THERMOELASTIC  STABILITY 

M.  Slemrod  and  E.  F.  Infante 

1,  Introduction 

Elastic  stability  is  usually  discussed  from  strictly  mechanical  con¬ 
siderations.  Recently,  however,  attempts  have  been  made  to  analyze  the 
influence  of  the  usually  neglected  thermodynamic  properties  of  elastic 
materials.  More  specifically,  one  may  ask  what  effects  the  second  law  of 
thermodynamics  has  on  the  asymptotic  stability  of  equilibrium  of  thermo¬ 
elastic  materials. 

KOITER  [1]  has  studied  the  general  nonlinear  thermoelastic  problem, 
and  for  materials  with  internal  friction  he  obtains  asymptotic  stability  of 
the  equilibrium  solutions.  ERICKSEN  [2]  has  posed  the  question  as  to  the 
asymptotic  stability  of  the  equilibrium  solutions  of  elastic  materials  with¬ 
out  imposing  the  assumption  of  internal  friction.  DAFERMOS  [3]  answered  this 
question  to  some  degree  by  obtaining  a  description  of  the  states  that  the 
material  approaches  as  t  -»  ». 

This  same  question  is  studied  here  in  a  more  general  setting  than  was 
done  by  SLEMROD  [4],  but  in  the  same  spirit!  it  is  shown  that  the  results  of 
[3]  can  be  obtained  as  a  simple  application  of  an  invariance  principle  for 
abstract  dynamical  systems  [4,5]. 

2.  Mathematical  Preliminaries 

The  principal  analytical  tool  to  be  used  is  a  generalization  due  to 
HALE  [5]  for  abstract  dynamical  systems  of  the  well  known  invariance 


principle  of  LASALLE  [6]  for  ordinary  differential  equations.  The  following 
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brief  presentation  of  this  tool  emphasizes  notation  and  concepts  to  be  used 
in  studying  the  problem  of  thermoelastic  stability. 

Let  R+  =  [0,«>)  and  let  9b  be  a  Banach  space  with  norm  ||4>||  for 
<t>  e  9b.  Then, 


Definition  2,1.  u  is  a  dynamical  system  on  a  Banach  space  9b  if  u 

is  a  function  u!  R  x  such  that  u  is  continuous,  u(0, $)  =  <t>, 

u(t+T, <t>)  =  u(t,u(r, <t>))  for  all  t,T  ^  0  and  all  <1>  in  9b.  The  positive 

orbit  0+  (<!>)  through  <t>  in  9b  is  defined  as  0+  (<i>)  =  (J  u(t,  <t>).  A  point 

teCT 

\|r  in  9b  is  an  equilibrium  point  if  0  (tJ/)  =  ijr. 

This  set  of  definitions  simply  generalize  familiar  notions  from  the 
theory  of  differential  equations  to  dynamical  systems. 

Definition  2.2.  A  set  M  in  ^  is  a  positively  invariant  set  of 
the  dynamical  system  u  if  for  each  <t>  in  M,  0  (<t>)  c  M. 

Definition  2.3.  A  set  M  in  ^  is  an  invariant  set  of  the  dynamical 

system  u  if  for  each  <t>  in  M  there  exists  a  function  U(s,<l>),  U(0,  $)  =  <t> 

defined  and  in  M  for  s  e  (-«>,<»)  and  such  that  u(t,U(s, 4>))  =  U(t+s,<t>)  for 
+ 

all  t  e  R  . 

Definition  2.2  is  well  known.  The  second  definition  is  used  to  extend 
backward  in  time  those  solutions  of  the  dynamical  system  which  lie  in  an  in¬ 
variant  set.  It  is  clear  that  if  a  set  M  is  invariant  it  is  positively 
invariant  but  the  converse  is,  in  general,  false. 


-3- 


Definition  2.4,  If  u  is  a  dynamical  system  on  &  and  V  is  a 
continuous  scalar  functional  on  &  ,  define  the  functional 

V(D)  =  TSJ  i  [V(u(t,4>))  -  V(4>)]. 
t-0 


Then 


Definition  2.5.  VI  -»  R  is  said  to  he  a  Liapunov  functional  on  a 
set  G  in  ^  if  Y  is  continuous  on  G,  the  closure  of  G,  and  if 
V(<l>)  ^  0  for  <t>  in  G.  Furthermore,  denote  by  S  the  set  S  =  in 
G|v(<t>)  =  0}  and  let  M  be  the  largest  invariant  set  in  S  for  the  dynamical 
system  u.  With  these  definitions  it  is  then  possible  to  prove; 

Theorem  2.1  (HALE  [3])«  Let  u  be  a  dynamical  system  on  & .  If  V 
is  a  Liapunov  functional  on  G  and  a  positive  orbit  0+(<i>)  belongs  to  G  and 
is  in  a  compact  set  of  !3)  then  u(t,  <t>)  M  as  t  -»  <». 

It  is  self-evident  that  in  applications  to  the  problem  of  asymptotic  sta¬ 
bility  of  an  equilibrium  point  i|r  it  is  necessary  to  show  that  M  =  {\|r}.  More¬ 
over,  it  should  be  emphasized  that  the  usefulness  of  this  theorem  in  applications 
depends  on  the  very  relaxed  assumptions  imposed  on  the  Liapunov  functional  V 
and  its  derivative  V.  These  conditions  should  be  compared  to  the  much  stronger 
conditions  imposed  by  standard  asymptotic  stability  theorems  (see,  for  example, 
PARKS  [7]). 

3.  Constitutive  Equations  of  Linear  Thermoelasticity 

A  material  point  is  identified  by  x  =  (x-^,x,_,,x^)  in  its  state  of 

TECftHTIOAL  LIBBAT'? 
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equilibrium  (no  stresses,  constant  temperature  =  TQ) •  The  displacement  field 
at  some  time  t  following  an  initial  disturbance  at  time  t  =  0  is  given  by 
u(x,t)  and  the  temperature  deviation  by  T(x,t)j  p(x)  denotes  the  density 
at  x  in  the  equilibrium  state. 

Let  ft  be  an  open,  bounded,  connected  set  in  E^  which  is  properly 
regular  [8];  let  8ft  denote  the  boundary  of  ft.  The  constitutive  equations 
of  thermoelasticity  are  taken  then  in  the  form 


pai  ■  (3a) 

pCD*  + 

where  body  forces  and  heat  sources  have  been  excluded.  In  these  equations 

^idk.0  ”  Cjik£  “  '"'kilo*  mid  ”  Kid  ~  K<3i  and  ^idk,^  mid  and 

K.  .  are  assumed  to  be  smooth  functions  of  x.. 
id  - 

Let  now  t  >  0.  By  a  classical  solution  of  the  mixed  initial-boundary 
value  problem  in  ft  X  (0,to)  we  mean  a  pair  (u,T)  satisfying  equations 
(3.1)  and  (3.2)  together  with  the  boundary  conditions 


u  =  0  on  8ft  x  (0,to)  (closed  boundary),  (3.3) 

T  =  0  on  8ft  x  (0,tQ)  (constant  temperature);  (3-4) 


and  with  initial  conditions 

(u(x,°),£(x,°),  T(x,0) )  =  (u^x),  u^x),  To(x)),  (3.5) 

where  u  (x),  u  (x)  and  T  (x)  are  given  functions  on  ft. 

— o  —  7  — o  —  o  — 


4.  The  Thermoelastic  Problem  as  a  Dynamical  System 

In  this  section  we  show,  by  recalling  some  results  of  DAFERMOS  [3], 
that  the  generalized  solutions  of  the  mixed  initial  boundary  value  problem 
described  above  can  be  viewed  on  an  appropriate  Banach  space  as  a  dynamical 


system  (ZUBOV  [9]).  Once  this  is  done,  the  application  of  Theorem  2.1  per¬ 
mits  us  to  draw  immediate  conclusions  on  the  asymptotic  Behavior  of  the  solu¬ 


tions  of  our  problem. 


( k)  ( k) 

Consider  the  Sobolev  spaces  Wr>  '(fl)  and  (^) j  4  =  1>2, ... 


(see,  for  instance  SOBOLEV  [10,  11],  AGMON  [12]).  Assume  that 

ess  inf  p(x)  >  0,  ess  inf  C^(x)  >  0, 
ft  n 

~  cikh>  ci>0  constant' 


(4.1) 

(4.2) 


(a  reformulation  of  the  Clausius -Duhem  inequality^  LANDAU  and  LIESHITZ  [13]) 
and  for  all  v^  £  ^20^^ 

;ncuiaTi,3Tk,/i  *  c2Vi,3vi,jd^  °2  >  0  c“Kt“<l,'5> 

a  general  property  of  the  tensor  of  the  elastic  modulii  for  infinitesimal  elas¬ 
ticity  (TRUESDELL  and  NOLL  [14]). 

Define  now  the  spaces  H  (G)  ~  W2o\^)  x  L  (ft)  X  Lp(ft)  with  norm 

l<vi,'ri,R>lo  *  ;f![wiwi  *  cijkfvi,jTk,f  ♦  ana  H(fl)  =  X 

( p  \  /  -i  \  °  onto 

¥^;(fi)  X  '(ft).  Define  the  map  P:  H  (ft)  ->  H-,(ft)  sending  (v.  w.  R)  e 

<—U  C-U  O  X  1  j  1  j 

H  (ft)  onto  (u.  v.  T)  e  H,(ft)  C  H(ft)  where  (u.,T)  e  W^(fl)  X  w£j?(ft)  is 

O  X  £-0  d{j 

defined  by  the  solution  ofi"  the  system 

4Ciokiuk,i8i,j—  *  -4tPwiei  + 

4V,P,id-  ’  ;fl[pci)E  +  mi3rovi,3^D  di 

for  every  D  e  W 2q4^)  .  The  mapping  P  is  linear,  well  defined  on  HQ( 2) 


and  one  to  one.  Hence,  defining  Pm=P°p° ...°  p  let  Hm(ft)  denote  the  range 
of  the  map  P^.  It  is  clear  that  P^  exists  and  maps  Hm(G)  onto  HQ(ft). 
Let  i|/  e  Hm(G)  and  define  1 1|/|  m  =  |  p"1^  .  Then* 
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Lemma  4.1  (DAFERMDS  [31) .  Hm  is  a  Banach  space  with  norm  |  *|m. 
HQ(ft)  D  H(ft)  D  ...  D  Hjfl)  algebraically  and  topologically.  Furthermore, 
is  dense  in  H^(ft)  for  m  >  £,  and  the  imbedding  I:  Hm(ft)  -» H^(ft) 
is  compact. 

Let  us  now  define  appropriately  a  generalized  solution  of  our  prob¬ 
lem* 


Definition  4.1.  (u.  ft.  T)  will  be  called  a  generalized  solution 

i,  ~  ‘  “ 

of  (3.1)  -  (3.5)  on  ft  x  (0,t  )  if  for  all  smooth  test  functions  (v.  R)  with 

o  1, 

compact  support  on  ft  and  vanishing  on  (1X0. 

t 

/  °/n( (4-t  )[pft.v.  -  C.  ..  .u.  v.  .  +  m.  .Tv.  . 

Jo  ft  v  o,L  i  x  idk i  k,jg  x,  0  id  x,d 


.  + 


PC  C 

+  — -  TR  .+  m.  .u.  (ft]  +  pft.v.  +  p  —  TR  + 

r0  id  1,0  i  i  rQ 

+  m.  .u.  R  -  — /^(K.  R  .)  .T  dtjdxdt  = 
id  i,  d  r  ov  id  ,i,0  - 


PC. 


=  -t  /n[pft  v.l  ,  +  — 2rp  rI  +  m.  .u  r|  ,  ,Jdx 

o  ftL  o.  x1 1=0  r  o  1 1=0  id  o.  .  ’  t=0  — 


l  '  'o 

With  this  definition  it  follows  that: 


i,d 


(4.4) 


Theorem  4.1  (DAFERMOS  [3]).  Under  assumptions  (4.1)  -  (4.3)  the 
triple  (u.  ft.  T)  describes  a  dynamical  system  on  H  (ft),  m  =  0,1,2,..., 

x,i,  m 

where  (u.  ft.  T)  is  the  generalized  solution  to  the  equations  of  linear 
thermoelasticity  satisfying  equation  (4.4).  Furthermore,  for  t  in  (0,tQ) 

*  l(ui  ,V  ,To>l« 

o’  o’ 


(*.5) 


n(m) 


where  Tv  '(Xjt)  denotes  the  generalized  mth  derivative  in  time  of  T(x,t). 


3.  Stability  Analysis 


The  problem  of  thermoelastic  stability  has  now  been  put  in  a  setting 
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appropriate  for  the  application  of  Theorem  2.1,  which  allows  us  to  obtain 

stability  results  in  a  simple  and  direct  manner. 

For  this  purpose,  fix  mil.  Then,  by  Theorem  4.1  and  (4.5)  it 

follows  that  for  any  initial  data  in  H  (ft)  the  trajectory  (u.  u.  T)(t) 

m  l,  x, 

will  lie  in  a  bounded  set  of  Hm(ft)  for  any  t  >  0.  Hence,  by  Lemma  4.1  the 
trajectory  remains  in  a  compact  set  G  of  Hg(ft),  i  <  m.  But  then  all  the  hy¬ 
pothesis  of  Theorem  2.1  are  met  with  =  Hg(ft).  For  simplicity  let  £,  =  0 

and  V  =  |  (u.  u.  T)|  .  From  (4.2)  and  (4.5)  it  immediately  follows  that 
x,  1,0 

•  1  2 

V  =  -  — fn K.  .T  .T  ,.dx  ^  -gJItIL  ,  C-,  >0:  therefore  the  set  S  is  the  set 

r  ft  id  ,i  ,3  -  3  l2  3  ’ 

S  =  {(u.  u.  T)  e  H  (ft)  I  ||t|L  =  0).  Let  now  M+  be  the  largest  positively 

i,  i,  o  L2 

invariant  set  in  S.  By  the  definition  of  generalized  solution  (4.4)  it 
follows  that  M+  =  {(u.  u.  T)  e  H  (ft)  |  ||t(  *,t)j|  =  0  for  t  1  0} .  Choosing 

1  y  1  p  O 

now  v.  =  0  in  (^.U)  it  follows  that  for  (u.  u.  T)  e  M+  it  is  necessary 


I  °fn  4rC(t-t  )R]m.  .u.  .dxdt  =  -t  /  m.  .u  R  dx: 

J  o  dt1  x  o'  J  13  1,3  -  oJn  13  o.  .  - 

3 

for  every  test  function  R.  Choosing  this  function  as  R(x,t)  = 

t-t 

o 

where  t](x)  is  an  arbitrary  test  function  and  u{ t)  is  the  C  tTbumpn 
function  of  Serrin  [15],  it  follows  that 

tis° 

for  (u.  u.  T)  in  M+.  Hence,  Theorem  2.1  applied  to  this  context  yields  the 

x,  1, 

desired  result* 

Theorem  5. It  For  any  initial  data  (u  ,u  ,T)  in  H  ,  m  g  1.  and 

l  l 

under  assumptions  (4.1)  -  (4.5)  (u^,tu,T)(t)  approaches  asymptotically  the  set 

((vi,vi,R)  e  wg)(ft)  x  L2(ft)  x  L2(ft)|  /nmiJ(x)vi  ..(x,t)dx.  =  fn±  .(x)u  (x)dx, 

9  J  °i^3 

t  S  0,  R  =  0)  in  the  norm  of  the  space  x  L2(ft)  x  L2(ft). 
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